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Abstract The Sun is a non-equilibrium dissipative system subjected to an en- 
ergy flow which originates in its core. Convective overshooting motions create 
temperature and velocity structures which show a temporal and spatial evo- 
3 ■ lution. As a result, photospheric structures are generally considered to be the 

direct manifestation of convective plasma motions. The plasma flows on the 
photosphere govern the motion of single magnetic elements. These elements 
are arranged in typical patterns which are observed as a variety of multiscale 
magnetic patterns. 

High resolution magnctograms of quiet solar surface revealed the presence of 
magnetic underdense regions in the solar photosphere, commonly called voids, 
which may be considered a signature of the underlying convective structure. 
■ The analysis of such patterns paves the way for the investigation of all turbulent 

convective scales from granular to global. 

In order to address the question of magnetic structures driven by turbulent con- 
vection at granular and mesogranular scales we used a voids detection method. 
The computed voids distribution shows an exponential behavior at scales be- 
tween 2 and 10 Mm and the absence of features at 5-10 Mm mesogranular 
scales. The absence of preferred scales of organization in the 2-10 Mm range 
supports the multiscale nature of flows on the solar surface and the absence of 
a mesogranular convective scale. 

Keywords: granulation mcsogranulation Magnetic fields, Photosphere 
1. Introduction 

The Sun is a non-equilibrium dissipative system subjected to an energy flow 
which originates in its core. A turbulent convective envelope is generated by this 
energy flow in response to the significant opacity of the outer solar region. The 
photosphere shows a state of reduced spatial symmetry due to the convective 
overshooting motions on the surface. The result of this reduced symmetry is the 
formation of temperature and velocity structures which show a temporal and 
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spatial multiscale evolution. As a result, photospheric structures are generally 
considered to be the direct manifestation of convective motions. 
Global plasma flows (i.e., differential rotation, meridional circulation and tor- 
sional oscillation) and the solar magnetic field, interacting with these turbulent 
convective motions, increase the complexity of emerging structures and of the 
star as a whole. The dynamics of the solar surface and its interaction with the 
magnetic field ultimately control the structure of the outer solar atmosphere and 
the hcliosphcre beyond. 

The plasma flows on the photosphere govern the motion of single magnetic 
elements. These elements are arranged in typical patterns which are observed as 
a variety of multiscale magnetic features, e.g., filamentary clusters or clumps of 
magnetic elements, magnetic network, clusters of facular points, active regions. 
The analysis of these multiscale magnetic patterns on the solar surface paves 
the way for the investigation of all turbulent convective scales from granular to 
global. Ycllc s Chaouche et al.l [2UTT1 report, indeed, that the magnetic field can 
provide a direct avenue to explore convective patterns since the magnetic flux can 
be measured directly using Stokes polarimetry techniques. 

Commonly, four different convective scales are identified on the solar surface: 
the granulation, which shows a typical horizontal scale of about f Mm and 
a short (a few minutes) lifetime, the mesogranulation, which shows horizontal 
length scales ranging from 5 Mm to 10 Mm and a lifetime of some hours, the 
supcrgranulation, which shows typical horizontal length scales from 20 Mm to 
50 Mm and a lifetime of about one day, and giant cells which show horizontal 
length scales of the order of lOOMm or larger. 

The wide range of values which are reported in the observed horizontal length 
scales of mesogranulation and supergranulation are most likely due to different 
methods of data analysis and choices of spatial and time filters. Even more impor- 
tant, however, is that this division, as suggested in Nordlund, Stein fc Asplund) 



2009, is largely of historical origin, and current evidence indicates that there is 
a continuous spectrum of motions, on all scales from global to sub-granular. 

In the last years, high resolution magnetograms of quiet solar surface (e.g., 
Berger et all IT9981 ISanchez Almeida! [20031 ILites et al.l 12008)) revealed the pres- 



ence of multiscale magnetic underdense regions in the solar photosphere, com- 
monly called voids, which may be considered a signature of the underlying 
convective structure. 

From magneto-convection simulations it results that at the granular scales the 
magnetic field elements are concentrated at the boundaries of granular cells (i.e., 
intergranular lanes). Several observations, i.e., |Dommguez Cerdena et al.|l2003"al 



IKhomenko et al.l 12003, Manso S ainz, Martinez Go nzalez & Asensio Ra mos) 120 
confirmed such numerical results. The presence of magneto-convective concen- 
tration at mesogranular scales has been predicted by numerical simulations 
HCattaneo et al.lBOOip and observed with 2D spectro-polarimetry ( jDommguez Cerdena et al| 



I2003H rYelles Chaouche et al.ll2UTTj), IR spec tropolarimetry ( [Trujillo BuenollM^ I 
and from magnetograms ([ Sanchez Almeidal I2003|) 



In more detail, Dommguez Cerdena et al ( 2003b assessed the existence of a web- 



like pattern with the spatial scale of mesogranulation by using a time sequence 
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of Inter-Network magnetograms observed at the solar disk center. The authors 
reported that the observed persistent pattern resembled a network with a spatial 
scale between 5 and 10 arcsec, which they identified as mesogranulation. 

More recently, Yclles Chao uche et aUf2011l confirmed the preferential location 
of magnetic elements in mesogranular cells using high spatial and temporal reso- 
lution IMaX-SUNRISE dataset. The mesogranular cells were identified using La- 
grange tracers driven by horizontal velocity fields which were computed via Local 
Correlation Tracking technique ([November fc Simon! I1988[) . In addition, the sta- 
tistical multiscalc analysis of the photospheric field geometry outlined the role of 
solar surface convective motions in the clustering and intermittent organization 



of magnetic flux elements in the quiet Sun (Lawrence, Ruzmaikin & Cadavid 
11993 | Uritsky fc Davila[ l2UTTj) . Moreover, |de Wijn fc Muller[ [20091 reported that 



it seems highly likely that the cells described by De Wijn et al.[ 120051 and the 



voids found bv lLites et al.[ |2"0"0"51 correspond to mesogranules. 

In order to address the question of which organization scales appear in mag- 
netic patterns formed by turbulent convective motions and how these scales could 
be revealed in high resolution magnetograms, we used a void detection method. 
Void detection methods are largely used in cosmology to study the voids from 
galaxy and clusters (jAikio fc M ahoncn, 1998]). These numerical methods define 
void structures of different particle distributions. In these methods galaxies are 
considered single particles. High-density particle regions surrounding low-density 
regions actually form the edges of these empty volumes, called voids. In the case 
of magnetic patterns on the solar surface we observe areas of low or absent 
magnetic activity surrounded by intense magnetic structures. Like the voids in 
the distribution of galaxies we define this underdense magnetic field regions in 
solar magnetograms as voids. 

In this work we employ a suitable automated void detection procedure to sin- 
gle out these underdense regions and to investigate their geometrical properties. 
We concentrate our analysis on the remarkable scales of granulation, the only 
length scale for which the convective signature is directly visible on the solar 
surface, and mesogranulation, whose physical origin is still under debate. See 
Nordlund, Stein fc Asplund , 2009 for a recent review on solar surface convection. 



In more detail, we are concerned with the study of distribution of voids in a quiet 
Sun high resolution magnetogram observed at disk center. 



2. Dataset 

We analyzed a high resolution magnetogram of a 302 x 162 arcsec 2 portion 
of the solar photosphere observed at disk center on 2007 March 10 between 
11:37 and 14:34 UT (FigO} The spectropolarimetric measurements were taken 
by the Solar Optical Telescope SOT/SP instrument aboard HINODE (?,?). The 
magnetogram has been derived from the Fel 630.15 nm Stokes I and V signals 
by adopting the center-of-gravity method (jRees Ik Semell ri979[) . 
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Figure 1. Hinode SOT/SP linc-of-sight magnetogram of a 302 X 162 arcsec 2 portion of the 
solar photosphere observed at disk center. 

We adopted the plate scale reported in the FITS header information of the Hin- 
ode SP, i.e., 0.1476arcsec/px and 0.1585arcsec/px along x and y, respectively. 
The (diffraction-limited) spatial resolution is about 0.32 arcsec. See lLites et al.l 
120081 for a detailed description of the dataset. 

3. Void detection algorithm and results 

In this work we adopted a 2-D version of the algorithm developed bv lAikio fc Mahonen[ [ 
119981 to detect voids in three dimensional red-shift galaxy surveys. Below, we 
formally define the void by analogy with |Einasto, Einasto fc Gramann| [1989 
and lAikio fc Mahonenl HM^l 

Let us consider a distribution of particles in a square region L 2 G 5i 2 with 
side length L. We define a scalar field, called distance field (DF), D : L 2 — > R 
as the distance of a given point x in L 2 to the nearest particle: 

D{x) =min{| x-X n |} (1) 

n 

where X n , n=l,. . ., N are the positions of the particles. 

The local maxima of the distance field, which are the points with the longest 
distance to the nearest particle, are the "centers" of empty regions, i.e., voids. 
To compute the DF we superimpose on the region L 2 a two-dimensional grid 
of k 2 elementary cells of equal area, where k = L/s. The parameter s, called 
resolution parameter, sets the spatial resolution of the void analysis. For each 
cell center we compute the distance to the nearest particle, thus obtaining the 
discrete DF D(x). 

Cells corresponding to the ith local maximum of D(x) arc identified with Mi. 
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Subvoid vm S L 2 is a region around Mi containing all the points x continuously 
connected to the maximum Mi along a monotonically increasing VD(i) path. 
More than one maximum of D(x) could fall inside a void; such local maxima 
belong to subvoids of the same void. Although a subvoid clearly belongs to a 
determined void, the subvoid does not necessarily cover the whole void area. In 
general we say that subvoids VMi and vm } belong to the same void V if: 

1. the maxima Mi and Mj are nearer to each other than the nearest particle 

di,j <mm{D{Mi),D{Mj)} (2) 

where d+j is the distance between Mi and Mj maxima, and D(Mk) is the 
value of the D(x) at point M k and ; 

2. if there exists a chain of subvoids {vM k } k=1 such that: 

d iA < min {D(Mi),D(Mi)}, d k , k +i < mm(D(M k ),D(M k+1 )} 
for all k=l,. . . , n-1 and 

d n ,j < min {D(M n ), D(Mj)} (i.e. there exists a chain of subvoids obeying the 
condition in EqJ2]). 
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Figure 2. Hinodc SOT/SP linc-of-sight two-levels magnctogram. The two-levels version of 
the magnctogram is computed taking into account only those pixels with magnetic absolute 
signal above the typical equipartition field strength of 200 G. 



The output of the algorithm contains the number N c of elementary cells in each 
void and the main radius Ry of the voids. The number of elementary cells defines 
the area of the void, in fact the area A is equal to N c multiplied the constant 
area s 2 of elementary cell, i.e., a pixel in our computational volume. The void 
length scale is computed as the equivalent diameter of a circular region with the 
same area, i.e., 2 x y/ A/ir, where A is the area of the detected void. 
The main radius Ry is the maximum of the D(x) inside the void, i.e., it is the 
maximum of the D(x) maxima D(M&) inside the void: 

R v = max W(M k )\ (3) 

M k £V 
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To start the algorithm we used a global threshold to get the two-levels version 
of the HINODE magnctogram. Since the mechanism of flux expulsion, for pho- 
tospheric convection, can concentrate magnetic flux up to the equipartition field 
strength, which is about 200 G at the solar surface ([Thomas! [19 90, Sola nki et al.l 
f 996), we decided to set the threshold value to 200 Gauss (absolute value). There- 
fore, in order to compute the two-levels version of the magnctogram we took into 
account only those pixels with magnetic signal above the typical equipartition 
field strength of 200 G. In FigJ3] we show the two-levels magnctogram. After 
detecting the voids present in the magnetogram (Fig[3]) we use an automated 
procedure to calculate area and position of each cell. The algorithm extracted 
1951 voids. In Figj5]we show the distribution of void length scale. 





Figure 3. Hinode SOT/SP line-of-sight segmented magnetogram. The image shows 1951 
identified voids. 



We applied the algorithm to several test images in order to validate its reliabil- 
ity and to study void distributions under controlled conditions. Particularly, we 
applied the void searching procedure on two ensembles of 20 simulated patterns 
of pseudorandom structures. The size of the computational volume and the 
number of structures were set equal to the size of the Hinode magnetogram and 
to the number of observed magnetic features, respectively The samples of the 
first ensemble did not include privileged size voids, while the samples forming 
the second ensemble included 100 pseudorandom round empty regions which 
generate voids of mesogranular scale. 

Final void distributions, for both ensembles, were computed by considering all 
the voids of 20 samples (Fig (5]). 

The final distribution shows two populations only for pseudorandom patterns 
of structures with non-overlapping mesogranular voids. The first population is 
identified by the 3 Mm maximum related to the correlation distance among 
structures. The second population is identified by the second maximum and it 
is due to the presence of the included mesogranular voids. 
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4. Discussion and conclusions 

As reported in paragraph [3J we applied the algorithm to a large SOT-HINODE 
magnetogram and the result is that 1951 magnetic underdense regions, i.e., voids, 
were identified. 

The calculated distribution of these 1951 voids (Figj4]) shows an exponential 
decrease between 2 and 10 Mm with a decay constant equal to 2.2 ± 0.2 Mm 
and a coefficient of determination R 2 = 0.98. The magnitude of R 2 indicates a 
very high degree of correlation. The absence of features at scales of 5-10 Mm 
indicates the lack of an intrinsic mesogranular scale. 

The effect of different threshold values is investigated from a comparison of corre- 
sponding void distributions. Our findings show that the decay constant decreases 
when the threshold increases from 120 Gauss to 260 Gauss. The distributions 
remain exponential with the same high degree of correlation, even though the 
mesogranular signature is not observed. 

Thus, we support the results reported bv lYelles Chaouche et al.ifeOlll that meso- 
granulation is not among the primary energy-injection scales of solar convection. 
In addition, these authors report that the Probability Density Function of the 
distance between magnetic footpoints shows a constancy of the slope at scales 
between 1 and 10 Mm and that the characteristic decay distance is approximately 
1.7 Mm. Moreover, they reported that only after about 20 minutes of integration 
a sharp mesogranular network appeared. 

Similarly, IBerrilli et al.l 120051 found that mesogranular photospheric features 
appear on the solar surface with a characteristic growth time of about 10 min- 
utes. They used the hexagonal normalized information entropy measure, H' hex , 
to observe the mesogranular scale emergence and to compute the associated 
time of formation. They found that the H' hex signal shows a maximum on the 
mesogranular scale of 8 Mm after 15-20 minutes of integration. 
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Figure 4. Void frequency distribution of 1951 identified voids. 
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Figure 5. Average void distributions for pseudorandom patterns of structures (dot) and 
pseudorandom patterns of structures with 100 round empty regions which generate voids of 
mesogranular scale (heavy continuum). In the second case, the final distribution shows two 
populations. The first population is identified by the 3 Mm maximum related to the correlation 
distance among structures. The second population is identified by the second maximum related 
to the mesogranular scale voids. 



The absence of specific scales of organization in the 2-10 Mm range supports 
the multiscale nature of flows on the solar surface, but does not rule out the 
possibility of convection structures with these dimensions. For example, large 
scale flows can be obtained as the result of the local merge of downward plumes 
in advective-interaction models (IRastl I2003|) . 

Exponential functions describe other physical quantities in turbulent convective 
systems, e.g., the probability density of normalized temperature fluctuations 
shows an interesting feature at the transition from soft to hard turbulent con- 
vection (jYakhotl I1989|) . As stated by authors, in the soft turbulent convective 
regime, when the Rayleigh number Ra < 10 7 , the probability distribution mea- 
sured at the center of the cell is Gaussian, while in the hard-turbulence regime 
(Ra > 10 8 ) the probability distribution is very close to exponential. The value 
of the Rayleigh number characterizing intensity of the turbulent convection in 
photosphere is Ra ~ 10 11 ( [Bray et al.| 11984)) which is consistent with the hard- 
turbulence regime. 

The hypothesis that there is a continuous spectrum of flows on the solar 
surface, from global to sub-granular ( Nordlund, Stein fc Asplund[ I2009[) . still 



needs to be thoroughly proven. With this in mind, we are developing a novel 
fast void detection algorithm to analyze a large dataset of MDI magnetograms 
acquired during the exceptional solar minimum at the end of solar cycle 23. 
The analysis of this dataset will allow us to investigate voids in the range of 
5-100 Mm in order to validate the above hypothesis. 
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